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Abstract 

The concern of this paper is the Cauchy problem for the Prandtl equation. This 
problem is known to be well-posed for analytic data [T^ITU]. or for data with monotonicity 
properties [HJ US]. We prove here that it is linearly ill-posed in Sobolev type spaces. 
The key of the analysis is the construction, at high tangential frequencies, of unstable 
quasimodes for the linearization around solutions with non-degenerate critical points. 
Interestingly, the strong instability is due to vicosity, which is coherent with well-posedness 
results obtained for the inviscid version of the equation [S]. A numerical study of this 
instability is also provided. 

1 Introduction 

One challenging open problem of fluid dynamics is to understand the inviscid limit of the 
Navier-Stokes equations 

' dtu" + ■ Vu" + Vp" - uAu" = 0, X G 17, 

V • u*" = 0, X G 17, (1.1) 

u^ldn = 0, 

in a domain 17 with boundaries, endowed with a no-slip boundary condition. Mathematically, 
the main difficulty is the lack of uniform bounds on the vorticity field, as the viscosity u goes 
to zero. In terms of fluid dynamics, this corresponds to a boundary layer phenomenon near 

dn. 

A natural approach to describe this boundary layer is to look for a double-scale asymp- 
totics, with a parabolic scaling in the normal direction. Consider the case 17 C M^. At least 
locally, any point x in a neighborhood of 317 has a unique decomposition 

X = yn{x) + x(x), x G 917, 

where y > 0, x is an arc length parametrization of the boundary, and n is the inward unit 
normal vector at 917. The velocity field can be written 

u''{t,x) = u''{t,x,y)T{x) + v''{t,x,y)n{x), 

where (r, n) is the Frenet frame. It is then natural to consider an approximation of the type: 

u''{t,x,y) PS u°{t,x,y) + u'^^{t,x,y/^/v), ^^^^ 
v''{t,x,y) « v^{t,x,y) + ^/l^v^^{t,x,y/^), 
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where 

u°(t,x,y) = u'^{t,x,y)T{x) + v'^ {t , x , y) n{x) 
satisfies the Euler equation with the no penetration condition, and 

describes a boundary layer corrector with typical scale in the normal direction. It is 
slightly more convenient to introduce 

u{t,x,Y) := u°{t,x,0) + u^^{t,x,Y), v{t,x,Y) := YdyV°{t,x,0) + v^^{t,x,Y). 

Indeed, inserting this Ansatz in the Navier-Stokes equations, we get formally 

' dtu + udxU + vdyu — dyu = [dtvP + u^dxU^) |j/=o, 1^ > 0, 

+ = 0, y > 0, 
(n,v) = o, y = o, ^^-^^ 

lim u = n'^L=o- 

These are the so-called Prandtl equations, derived by Ludwig Prandtl [12]. Note that the 
curvature of the boundary does not appear explicitly in the system. It is however involved 
in ()1.3p through the Euler field, and through the interval of definition of the arc length 
parametrization x. Up to our knowledge, all studies deal with one of the three following 
cases: x S M, x G T, or < x < L, supplemented with a condition on it at x = 0. The first 
and second choices are convenient to describe phenomena that are local in x. The case x G T 
may also model the outside of a bounded convex obstacle. Finally, the third configuration is 
adapted to the spreading of a fiow around a thin obstacle, where x = corresponds to the 
tip of the obstacle. 

Although this formal asymptotics is very natural, its validity is not clear. As emphasized 
by physicists, including Prantdl himself, it may not hold uniformly in space and time. One 
reason is the so-called boundary layer separation^ which is observed for fiows around obstacles, 
see [^. Nevertheless, the description (II. 2p fits with many experiments, upstream from the 
separation zone. In any case, to understand the relevance and limitations of the Prandtl 
model is a crucial issue. 

From the mathematical point of view, one must address two problems: 

1. The well-posedness of the Prandtl equation. 

2. The justification of the expansion (II. 2p . 



These two problems depend crucially on the choice of the underlying functional spaces, es- 
pecially on the regularity that is required in the tangential variable x. Indeed, the main 
mathematical difficulty is the lack of control of the x derivatives. For example, v is recovered 
in (|1.3p through the divergence condition, and in terms of x-regularity, behaves broadly like 
dxU. This loss of one derivative is not balanced by any horizontal diffusion term, so that 
standard energy estimates do not apply. 

Within spaces of functions that are analytic in x G M, 1" G Sammartino and Cafiish 
have overcome these problems, justifying locally in time the boundary layer asymptotics 
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[13^ I14j . But for more "realistic" functional settings, the way solutions of (jl.ip behave is 
still poorly understood. Various instability mechanisms, that are filtered out in an analytic 
framework, become a huge source of trouble. For example, when the viscosity is small, the 
Navier-Stokes equation admits exponentially growing solutions which are both small-scale 
and isotropic in x, y. Their evolution is lost in the anisotropic Prandtl description. 

This remark was used by Emmanuel Grenier in [6j, who relied on the so-called Rayleigh 
instability for inviscid flows to show that the asymptotics (jl.2p does not (always) hold in the 
Sobolev space (see [6] for a precise statement). However, the relevance of this asymptotics 
in LP, or its relevance in the absence of Rayleigh instabilities, are still open issues. 

Above all, the local in time well-posedness of the Prandtl equation for smooth (say Sobolev) 
initial data has been so far an open problem. Up to our knowledge, the Cauchy problem has 
only been solved in two settings: 

i) X S M, with data that are analytic in x : see \13\ [TO] for more. 

ii) < X < L, with data that are monotonic in y : see [TS] for more. 

One may also cite article in which blow up in time of some smooth solutions is exhibited. 
Finally, let us mention the interesting work [8], in which the inviscid version of (II. 3p is 
analyzed (no dyU in the equation). Interestingly, for a smooth initial data, this equation 
turns out to have an explicit solution through the method of characteristics. In particular, 
starting from a smooth data, one recovers locally in time a smooth data. More precisely, 
there is only a finite loss of x-derivatives, so that the Cauchy problem is (weakly) well-posed. 
We refer to |8j for all details. See also papers O [1] on the hydrostatic equations, that share 
some features with Prandtl equations. For more on Prandtl equations, see the review [3]. 

On the basis of the inviscid result, it seems reasonable to bet for well-posedness of the 
Prandtl equation ()1.3p in Sobolev type spaces. The aim of this paper is to show that it is 
actually linearly ill-posed in this framework. As we shall see later on, the reason for ill- 
posedness is a strong destabilization mechanism due to two ingredients: viscosity, and critical 
points in the base velocity profile. In particular, it does not contradict the positive results 
obtained in the inviscid case and for monotonic data. 

Let us now describe our results. We restrict our selves to {x,Y) e T x M+, and u° = 0. 
To lighten notations, we write y instead of Y. The Prandtl equation comes down to 

' dtu + udxU + vdyU — dyU = 0, in T x M"*". 

dxu + dyv = 0, inTxM+, (1.4) 
('u,'t;)L=o = (0,0), lim n = 0. 

Let Us = Us{t,y) a smooth solution of the heat equation 

dtUs-dyUs = 0, Us\y=o = 0, Us\t=o = Us, (1.5) 

with good decay as y ^ +00. Clearly, the shear velocity profile {us,Vs) = {us{t, y),0) satisfies 
the system (|1.4p . We consider the linearization around {us,Vs), that is 

' dfU + UsdxU + vdyUg — dyU = 0, in T x M+. 

, dxU + dyV = 0, in T X M+, (1.6) 

('u,'t;)L=o = (0,0), lim u = 0. 
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We wish to study well-posedness properties of (jl.6p . for a certain class of velocities Ug. In 
this view, we introduce the following functional spaces: 

W/^--(M+) := {/ = /(y), e-VGt^^'°°(M+)}, Va, s > 0, 
with 11/11 vi/='°° := \\e°'^f\\w''-°^, and 

Ea,0 := |^x = n(x,y) = J^n'=(y)e*^^ ll^i'll^'- < C«,/3 e'^^l'^l, VA; I , Va,/5>0, 

with ||n|b^_, := sup,e/5W||^'=||^o,oo. 

Note that the functions of Ea,i3 have analytic regularity in x. They have only L°° regularity 
in y, with an exponential weight. More regularity in y could be considered as well. Let 
a, /? > 0. We prove in the appendix the following result: 

Proposition 1 (Well-posedness in the analytic setting) 

Let Us G C°(M+; Wi'°°(M+)). There exists p> such that: for all T with f3 - pT > 0, 
and all uq G Ea^js, the linear equation (11.60 has a unique solution 

U £ C {[0,T); Ea,l3~pT) , u{t,-) £ Ea,l3-pt, u\t=o=uo- 



In short, the Cauchy problem for ()1.6p is locally well-posed in the analytic setting. We shall 
denote 

T{t,s)uo := u{t,-) 

where u is the solution of (|1.6p with u\t=s = uq. As the spaces Ea^p are dense in the spaces 

:= /7™(T^,TyO'°°(M+)), m > 0, 
this makes sense to introduce the following notation: for all T G C{Ea^f3, Ea^p'), 

\\Tuo\\h'^2 



\'E\\c(H"^i,H"'2) — sup 



UQ 



that belongs to M-(_ U {+cxd}. In particular, it is infinite when T does not extend to a bounded 
operator from H"^^ to H"^^ . The main result of our paper is 

Theorem 1 (Ill-posedness in the Sobolev setting) 

i) Let Us e C°(E+; Ty^'°°(M+)) n (M+; Ty^'°°(M+)). Assume that the initial velocity has a 

non- degenerate critical point overM'^. Then, there exists o" > 0, such that for all 6 > 0, 

sup \\e-''^^~''^'y^\ T{t, s)\\c(H--,H---^) = +00, Vm > 0, /i G [0,1/2). 

0<s<t<5 

ii) Moreover, one can find solutions Us of (II. 5p and a > such that: for all 5 > 0, 

sup ||e-'"(*-")\/^r(t,s)||^(j^™i_^„.2) = +CX), Vmi,m2>0. 

0<s<t<S 
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This theorem expresses strong hnear ih-posedness of the Prandtl equation in the Sobolev 
framework. It is a consequence of an instabihty process, that holds at high tangential fre- 
quencies. We will show that some perturbations with tangential frequency A; S> 1 grow like 

The outline of the paper is as follows. Section [2] gives a formal description of the instability 
mechanism. It relies on an asymptotic analysis of ()1.6p . in the high tangential frequency 
limit. Thanks to this analysis, we show that ill-posedness for the PDE ()1.6p comes down to 
a "spectral condition" for a reduced ODE, namely: 

(SC) There exists r G C with Xmr < 0, and a solution W = W{z) of 

{r-zr^W +^^{{r-z')W)=0, (1.7) 
such that lim W = 0, lim W = 1. 

2— > — OO 2— > + 00 

This spectral condition is studied in section [3l and shown to be satisfied. On these grounds, 
we prove theorem [H cf section HI We end up the paper with numerical computations, which 
emphasize that our instability mechanism is effective. 



2 The instability mechanism 

In this section, we describe the destabilization of system (II. 6p . leading to the ill-posedness 
theorem. As we shall see, it takes place at high tangential frequencies, say 0(l/e), and has a 
typical time 0{^/e). At this timescale, the time dependence of the base velocity {us{t,y),0) 
will not play an important role. Thus, to understand the instability mechanism, we can 
consider the simpler equation 

' dtu + Usd^u + vU', - dlu = 0, in T x M+. 

< d^u + dyv = {), inTxM+, (2.1) 

{u,v)\y=o = {m- 

Handling of the real equation, that is with Ug instead of C/^, will require minor modifications, 
to be achieved in section HI 

System (|1.6p has constant coefficients in t and y, so that we can perform a Fourier analysis: 
we look for solutions in the form 

u(t,y) = e^'=('^('=)*+^)n^'(y), = A;e*'=(^(^')*+^){;'=(y), A; > 0. (2.2) 

As we are interested in high tangential frequencies, we denote e := 1/k <^ 1, and write 
uj{e) instead of u}{k), Us{y),V£{y) instead of (y) , v'' (y) . The divergence condition yields 
v'^{y) = —ius{y). Using this relation in the first equation in p.ip . one ends up with 

{u{e) + Us) < - Uive + i = 0, y > 0, 

V£\y=0 = V'^\y=0 = 0. 

Thus, the high frequency limit e — > in variable x yields a singular perturbation problem in 
variable y. To investigate this problem, one must first consider the inviscid case e = 0. 
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2.1 The inviscid case 

When e = 0, one can a priori only retain the impermeabihty condition. The appropriate 
problem is 

{oj + Us) V - = 0, y> 0, v\y=Q = 0. (2.4) 

This spectral problem, as well as the corresponding evolution equation, have been studied 
exhaustively in [8]. Clearly, there are non-trivial solutions if and only if oj belongs to the 
range oi —Ug- Moreover, the couples 

LOa = -Us{a), Va = H{y - a) (Us - Us{a)) , a > 

where H is the Heaviside function, satisfy (j2.4|) . Note that the regularity of Va depends on the 
choice of a. When a is a critical point, it belongs to Wa^{^'^) with a discontinuous second 
derivative. Otherwise, it is only in Tyi'°°(R"'"), with a discontinuous first derivative. Luckily 
enough, the additional boundary condition v'Jy=o = is also satisfied. 



2.2 The viscous perturbation 

When e is not 0, the inviscid eigenelements Ua, Va do not solve ()2.3p . All boundary conditions 
are satisfied, cf the above remark, but the equation is not. First, there is a 0(e) remaining 
term for y > a. More importantly, Va is not smooth at y = a, whereas a solution of this 
parabolic equation should be. 

Nevertheless, at least if a is a non-degenerate critical point, there is an approximate 
solution near {iOa,Va)- We shall establish this rigorously in section SI We just give here a 
formal expansion. It reads 

< 

Veiy) 

where r G C, and V = V{z) quickly tends to zero as z ^ ±00. Note that the approximation 
of Vg has two parts: the "regular" part 

vl'^iy) = H{y-a) {Us{y) - Usia) + e^l\) 

and the "shear layer part" 

For uj{e) = —Us{a) + e^/^r, the function vl^^ solves (12.30 up to 0(e), away from the critical 
point y = a. However, it has a jump at y = a, together with its second derivative. The role 
of the shear layer vf, which concentrates near y = a, is to cancel these discontinuities. Still 
formally, we obtain the system satisfied by the profile V: 

' (^r + f/^'(a)0y' - Ui'ia)zV + iV^-^^ =0, z / 0, 

limy = 0. 

^ ±00 



(2.5) 
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Let us point out that this system is a priori overdetermined, as jump and boundary conditions 
provide too many constraints. This justifies the introduction of the parameter r in the Ansatz 
(I2.5p . As we shall see below, there is a r for which system (12. 6p has a solution. Moreover, 
Tmr is negative. Hence, back to the Fourier representation (|2.2|) . the k-th mode will grow in 
time like e^*. This is the key of the instability mechanism. 

To see how the condition (SC) of the introduction steps in, we need a few rewritings. 
First, r + C/"(a)%- satisfies the equation in ()2.6p . We therefore introduce 



Viz) 



which leads to 



r + U'Jia)- 



V{z) + l^^[r + U'J{a)- 



V' - U'^{a)zV + =0, ze 



\\mV 

— oo 



T + U':{a)- 



0, V ~+oo 

Then, we introduce W such that 

V = {t + U'^{a)z^/2) W. 



We get: 



[r + U':{a)z^/2fj-W + ^^^[[r + U:{ 



a) z' 



72) W) = 0, 



\miW = 0, 



limH^ 

+ 00 



1. 



Finally, we perform the change of variables 



V2 



U"{a)\ 



1/2, 



2i/4|[/"( 



Dropping the tildes leaves us with the reduced ODE 

(r+ sign{U':{a))z'f^W + i-^{{r+ sign(t/^'(a)) z^) w) = 0, 
limW = 0, limTy = 1. 

— oo +00 

If U'J{a) < 0, it is exactly the system in (SC). If on the contrary U'J{a) > 0, and if {t,W) 
satisfies the system in (SC), then (r := — r, W := W) satisfies the above system. In both 
cases, back to the original system (12. 6p . condition (SC) gives a solution (r, V) with Imr < 0. 
In particular, this ^/£ correction to the eigenvalue is a source of strong instability, leading to 
ill-posedness. 

The proof of Theorem [H which is based on this formal shear layer phenomenon, is post- 
poned to section m In the next paragraph, we focus on condition (SC), and prove that it is 
satisfied. 
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3 The spectral condition (SC) 

We need to study the existence of heteroclinic orbits for the ODE (|1.7|) . Note that = 1 is 
a solution. Equation ()1.7p can be written as a second order equation in X = W': 

i{T-z^)X" - 6izX' + {{t - z'^f - 6i) X = 0. (3.1) 

To show that (SC) holds, we proceed in three steps. 

Step 1. We consider an auxiliary eigenvalue problem: 

For its study, we introduce the weighted spaces 

£2 := |n E L^, ^ {z^ + lf\u? < +cx)| , 

^1 := [u G I {z^ + l)Vl' + l{z^ + ^f\uf < +00 1 C\ 

I Jk Jm. J 



with their obvious Hilbert norms. We see A as an operator from D{A) := |u € Tl^, Au G 
into C^. Our goal is to show that A has a positive eigenvalue. 

By standard arguments, the domain D{A) is dense in C^. Moreover, for any u in D{A), 
there is a sequence n„ of smooth functions with compact support, such that u„ — > u in 
and Aun — > Au in Integration by parts and use of this density property give easily that 
A is symmetric, i.e. 

yu,veD{A), {Au\v)c2 = {Av\u)c2 (3.3) 
and that for A large enough 

i{X-A)u\u)c2 = X [ (z2 + l)4|n|2-6 [ (z2 + l)2|^|2^ f (^2^1)3|^'|2 > (3.4) 
JR JR JR 2 



Then, the coercivity condition (13. 4p allows to apply the Lax-Milgram lemma. It implies the 
invertibility of A — ^, with 

\\iX-Ar'f\\c2 < WiX- Ay' fWn. < C\\f\\c2. 

Moreover, from (|3.3|) . (A — A)~' is selfadjoint, and so is A. 

First, let us prove that A has positive spectrum. To do so, we claim that it is enough to 
find u G DiA) with [Au \ n)£2 > 0. Indeed, suppose a contrario that (7{A) is contained in 
R_. Then, by the spectral theorem, 

Vq > 0, ||(^ - a)-'\\ = - — 3—— < a-\ 



We deduce: for all u G D{A) 



\u\\j-2 < ||(^ - a)u\\l2 
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Expanding the scalar products, we obtain 



< a~^\\Au\\l2 - 2a~^{Au\u)c2 

In the Hmit a — > +00, we get {Au \ n)£2 < for all u £ D{A). This proves our claim. Prom 
there, we simply take u = e"^^ . A straightforward computation gives 

/ . , s 439 ^ 
{Au u) = —V^ > 0, 
biz 

and so a (A) has a positive subset. 

It remains to exhibit a positive eigenvalue inside this positive subset of the spectrum. We 
remark that the operator A can be split into 

.... 1 // 62; , . 6 

A = Ai + A2, Aiu := u + , p , ^.^ u , A2U := u. 

+ 1 (z^ + 1)^ (z^ + 1)^ 

On one hand, the operator Ai is negative, and by Lax-Milgram Lemma, Ai — A is invertible 
for any A > 0. Thus, cr{Ai) C M_. On the other hand, let Un a sequence with Un and AiUn 
bounded in C^. This implies that ii„ is bounded in Ti^, and so has a convergent subsequence in 
Lf^^. Moreover, is equi-integrable over M. Finally, it implies that A2Un has a convergent 
subsequence in C^, which means that A2 is Ai— compact. Hence, the essential spectra of A 
and Ai are the same, see In particular, the positive part of cr(A) is made of isolated 
eigenvalues with finite multiplicity. Eventually, we state: there exists a > 0, and u in D{A) 
satisfying (13. 2p . 



Step 2. We wish to convert the eigenelements (a, u) of the previous step into an appro- 
priate solution (r, X) of (13. We set f = — a^/^, and z = a~^^^z, Y{z) = u{z). Dropping 
the tildes, we obtain a solution of 

{t-z^)Y" - 6zY' + {{T-z^f-6)Y = 0. (3.5) 

By a classical bootstrap argument, Y is smooth. Moreover, it inherits from u its integrability 
properties at infinity. Actually, the behaviour of Y can be further specified, as shown in: 

Proposition 2 The function Y admits a unique extension, still denoted by Y , that is holo- 
morphic in z and satisfies (13. 5p in the simply connected domain 



Ur := C\(^ -ioo, -i|r|i/2 U i|r|i/^ +ioo ) 



Moreover, in the sectors arg z € (— 7r/4 + 5, 7r/4 — 5) and arg z G (37r/4 + 5, 57r/4 — S), 6 > 0, 
it satisfies the inequality 

\Y{z)\ < Cs exp(-zV4). 

Proof. This proposition follows from the general theory of ODE's with holomorphic coef- 
ficients. The existence of a holomorphic solution is well-known, because the coefficient r — z^ 
does not vanish on Ut- As regards the inequality, we rewrite equation p.ip as the first order 
system: 

ly = zA.)y. y={.-^iy). = (3.6) 
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In particular, A is holomorphic at infinity, with ^(oo) = ( 5 J )• distinct eigenvalues 

±1, with eigenvectors (±1) • 

Hence, we can apply [21 Theorem 5.1 pl63]: in any closed sector S inside which Rez^ does 
not cancel, there exists solutions y± (depending a priori on S) with the following asymptotic 
behaviour as \z\ — > +00: 

where a± is a complex constant, and P±{z) is a polynomial of degree 2. Moreover, the leading 

2 

term of P± is ib^. Following the same scheme of proof, we get in the present case: 

y± = (±1) , «± = -li±r + 7), P± = ±zV2. (3.8) 

As our solution Y is integrable over M, it is necessarily proportional to the decaying solution. 
The bounds in proposition [2] follow. 

Now, as Y is defined on Ur, we can perform the complex change of variable: 

z := e^"/^z, f := e^'^V 

Note that f has a negative imaginary part. Moreover, for z real, the original variable z belongs 
to the sectors arg z £ (— '/r/4 + 5, 7r/4 — 5) or arg z £ (37r/4 + 5, 57r/4 — 5), with 5 = 1/16. By 
the proposition [21 the function X{z) = Y{z) satisfies the estimate |X(-S)| < C exp(— l/4z^). 
Finally, dropping the tildes yields a solution r, X of (13. ip , where X decays at infinity. This 
concludes step 1. 

Step 3. To deduce from the previous step that (SC) holds, it is enough that J^X(z)dz be 
non zero. If so, one can consider 

-1 



W{z) := ( [ X{z')dz') r X{z')dz' 
\Jr J J-00 



which clearly satisfies all requirements. 

Let us assume a contrario that J^X{z)dz = 0. Then, the function 



Viz) := {t-z') / X{z')dz' 
J —00 

is a solution of 

{t-z^)V' + 2zV + iy(^) = 0, 

which decays exponentially as z goes to ±00, together with all its derivatives. Differentiation 
of the equation gives 

{t-z^)V" + 2V + iV^^^ = 0, 

Then, we multiply by V", that is the complex conjugate of V" , and integrate over R. Simple 
integrations by parts yield: 

(^r-z')\V"\^ - 2 / \V'\^-i [ = 0. 



10 



The imaginary part of this identity yields 

Imr [ \V"\'= f |y(3)|2 
Jr Jr 

which contradicts the fact that Imr < 0. Thus, the condition (SC) is satisfied. 

4 Proof of ill-posedness 

Theorem [T] win be deduced from the formal analysis of section [2l This analysis was performed 
on (|2.ip . in which possible time variations of u^were neglected. To account for the original 
system p.6p will require a few modifications, notably in the choice of the approximation ()2.5p . 
We will distinguish between the parts i) and ii) of the theorem. 

4.1 Ill-posedness for general Ug 

Let Us satisfying the assumptions of part i). Let a be the non-degenerate critical point of 
Us\t=o = Us- For the sake of brevity, we consider the case C/"s(a) < 0, the other one being 
strictly similar. The differential equation 

dtdyUs{t,a{t)) + d^Us{t,a{t))a'{t) = 0, a(0) =a 

defines for small time t < to a non-degenerate critical point a{t) of Us{t,-). Let then t,W 
given by condition (SC). We set 

V := {t- z^)W - 1m+ (r - z") . 

In the light of section [2l we introduce, for e > and t < to: 

u;{e,t) := -Us{t,a{t)) + -j^\d^Us{t,a{t))\^/^r 
as well as the "regular" velocity 

vl'S{t,y) := H{y-a{t)) (^Us{t,y) - Us{t,a{t)) + '-^\d'yUs{t,am'/\^ , 

and the shear layer velocity 

vl^{t,y) := '-^^{y - a{t)) |ggn,(t, a(t))| y (|agn,(t, a(t))| (^^^^^^^ ^ , 

where (p is a smooth truncation function near 0. We then consider the following velocity field: 

n,{t,x,y) := e-"^C/,(t, y), C/,(t,y) = i e^" /o -(M^^ a, {vT'{t,y) + vf{t,y)) , 
v,{t,x,y) := e^^-'^V,{t,y), = e"! e^"' /o -(-.^)'^^ (vl'%t,y) + vl\t,y)) . 

In order to have a field that is 27r-periodic in x and growing in time, we take e := ^, with 
n G N*. One verifies easily that Ug = e*^ ^^Us{t,y) is analytic in x, and VF^'°° in t,y. 
Moreover, we have the bounds 

ce^ < \\Ue{t,-)\\^2^oo < Ce^, (4.1) 
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for positive constants c, C and do that do not depend on e. 

Inserting the expression for , f e into the hnearized Prandtl equation (|1.6|) , we obtain 

' dtUe + UsdxUe + v^dyUs - dyUs = r^, in T x M+. 

dccUs + dyVe = 0, in T X M+, (4.2) 
iu,v)\y=o = (0,0). 

The remainder term reads = e*^ ^^Rs{t,y), with 



,2, 



i?,(t,y) = e--'/o-(^'^)'^^(-e-i(n,(t,y) - n,(t,a(t)) - a2n,(t, a(t))|-) 5,<'(t, y) 

+ e-^(dyUs{t,a{t)) - dlus{t,a{t))y)vf{t,y) 
- iedlvT^it^y) + idtdy(vT3{t,y) + <'(i,y)) + 0{e^) 

The 0(e°°) gathers terms with derivatives of (p: as the shear layer profile V decreases ex- 
ponentially, and the derivatives of (p{- — a) are supported away from a, their contribution is 
indeed exponentially small. Straightforwardly, 

\\Re{t,-)\\^o,^ < Ce^, (4.3) 

with the same ctq as in (j4.ip . 

We are now in a position to prove part i) of Theorem [H Let us assume a contrario that 
for all o" > 0, there exists m > 0, G [0, 1/2) and 5 > such that 



0<s<t<5 

Let 



sup ||e"''^*~'VI<^-lT(t,s)||£(H-,//— M) < +00. 

0<s<t<(5 



T,{t,s) : <'°°(M+) ^ <'°°(M+) 

the restriction of T{t,s) to the tangential Fourier mode e~^. Namely, T{t,s) ^e*^ C/q 

g«£ J|,(t^ s)[/q. Similarly, we denote Le = e~*^ ^^Le*^ where L is the linearized Prandtl 
operator around ti^. We have, for all < s < t < (5, 



ll^(*;'5)ll£(Vl/°'°°) 



o-(t-a) 



Let U = U{t, y) the solution of dtU + L^U = 0, that coincides initially with the approximation 
Us- On one hand, we get 

\\U{t,-)\\^o.o. < Ce-^e^\\U{0,-)\\^o,^ < C'e-'^e^. (4.4) 

On the other hand, the difference U = U — satisfies, for all t < S, 

U{t,-) = I Ts{t,s)Re{s)ds. 
Jo 
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Estimate (|4.3p implies 

\\U{t,-)\\^o,oo < Ce-f" / e^^ev^ds < C'e^/^-^gv^^ 

^0 

as soon as o" < do. Combining this with the estimate (|4.1|) . we obtain the lower bound 

\\U{t,-)\\^u,^ > \\Ue{t,-)\\y^o,^ - \\lJ{tr)\\wO,^ 

> ce - Ce'^-i/^gv^ 

For e small enough, we get 

\\U{t,-)\\^,,^ > c' e^, 

which contradicts the upper bound (j4.4|) . as soon as o" < ctq and t I V^- This 

achieves the proof of part i) . 

4.2 Stronger ill-posedness for specific Ug 

It remains to handle part ii) of Theorem [TJ Roughly, we must find some Ug for which 

e~"\f\^\^*~'^'> T{t,s) fails to be bounded from H™- to H"^~^^, /x > arbitrary. Using no- 
tations of the previous paragraph, the keypoint is to build, for any N, a growing solution 

Ue,N of 

dtUe,N + LeUe,N = Re,N, where \\Re,N{t,-)\\wO,^ < Cn [e^ + 1^^) e^. 

Indeed, we can then take N + 1/2 > /i, and conclude along the same lines as above. 

So far, we have not managed to improve the approximation of the previous paragraph for 
general Us- This explains the technical restriction fi G [0, 1/2) of part i). In order to obtain 
a refined approximation, we consider some special profiles: we assume that Us{0,y) = Us{y), 
for some exponentially decreasing Ug, satisfying in the neighborhood of a > 0: 

Us{y) = C/"(a)^^^, U':{a)<0. 

Notice that o is a non-degenerate critical point of Ug- For such profiles, the approximation of 
the previous paragraph reads 

Ueit,y) = ze^^-'/o-CM-^^aj, {vl'nt,y) + vf{t,y)) , e = i, n G N.. 

Using that Ug is quadratic near y = a, one can improve this approximation through an 
expansion of the type 

TV 
i=l 

with additional terms vl'^^^ . Let us briefly explain the construction of these extra terms. The 
error terms due to Ue divide into three categories: 
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1. Shear layer terms involving derivatives of As mentioned before, they are O and 
require no correction. 

2. Terms that come from the replacement of Ug by its Taylor expansion in the shear layer 
equation. They read 



Re,2 

We write 



'(dyUsit,ait)) - dlus{t,a{t))y^vf{t,y). 



/ 



y (z-a(t))\ ^ ^ , 

^ ^^^dlus{t,z)dz 

it) 2 



dyv'J{t,y) 



^ ry 

Ja{t) 



y {z \d^d^yUs{0,z)\ dz \dyvl^{t,y)\ + 0(t2^)e^ 



,zni fy {z-a{t)y 



ry 

Jait) 



k=o 

2N-1 



E It d'y+^'Usiz) dz dyvf{t,y) + 0(i2^)e^. 



k=0 



k\ 



The second inequality stems from a Taylor expansion of Ug with respect to t. As Ug 
satisfies the heat equation, each time derivative can be replaced by two space derivatives, 
so the third line. Because Ug is quadratic in a vicinity of a{t) (for short times), and vf 
and its derivatives decay exponentially fast, we end up with 

A similar bound holds for Re,2- Hence, these remainders do not require correction. 

3. Terms that come from the time derivative and the diffusion. We focus here on the time 
derivative, as the diffusion term is simpler and has smaller amplitude. This is 

Re,3 := ie''~' ^o-'i-^^)ds d^dy(vl''\t,y) + vf{t,y)). 

Proceeding as for -^£,2) that is with Taylor expansions in leads to 

As regards the regular part, 

gi£-i S'^oj{e,s)ds dtdyvl^>^{t, y) = e'^~' fo^(''')ds _ a(t)) dtdyUg 

= e'^~' fo'^^^'^)'^^ H{y - ait)) Fit,y) + O^^t^iVe^j 



where 



2N-1 

fe=o 
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comes again from a Taylor expansion in t. The nice thing about this 0(1) term is that it 
is identically zero in the vicinity of y = a{t) (for short times). As a result, the Heaviside 
function H{y — a{t)) in front of it does not create any discontinuity, and no extra shear 
layer term is necessary. One takes care of this source term by the introduction of 

v'/^^ = H{y-a{t))(us{t,y)-co{e,t)) T ^^^''\ dz, 

^ ' Ja(t) {Us{t,z) - Uj{£,t)) 

SO that 

[/I := -iee'^~' ^o^''''')'^'dyvl''^^3 

solves 

dtU^ + L,U^ = ie^^"'/o-(-'^)'^^/7(y-a(t))F(t,y) + 0(e). 

Proceeding recursively, we obtain an approximation as accurate as we want. This ends 
the proof of the theorem. 

5 Numerical study 

In this last section, we present numerical illustrations of the instability process. 
5.1 Numerical test of (SC) 

To check (SC) numerically, it is more convenient to reformulate it with an Evans function. 
We know from Step 2, section (Sj that there are solutions y±{z) of ()3.6|) satisfying (I3.7p - ()3.8p . 
Back to the ODE (jS.ip . this yields independent solutions X+(r, •) and X_(r, •) respectively 
growing and decaying as z goes to +oo. Furthermore, the following asymptotics holds: 

IT 7 1 

X±{t,z) ~ exp(±-Az2), 
d,X±iT,z) ~ ±Az±5x-iexp(±^Az2), 

with A = Thus, the functions 

W.{t,z) := / X_iT,s)ds, W+iT,z) := / X+{T,s)ds, Wo{t, z) = 1, 
Jz Jo 

seen as functions of z, form a basis of solutions of ()1.7p . They are respectively decaying, 
growing and constant at +oo. As ()1.7p is preserved by the change of variable z —z, the 
functions W-{t, —z), W+{t, —z), Wq form a basis as well. They are respectively decaying, 
growing and constant at —oo. The existence of the heteroclinic orbit is the same as the 
existence of some constants A and B such that 

1 = AW-{-z,t) + BW-{z,t) 

for all z, or equivalently 
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This last condition is easily seen to be equivalent to 

W-{t,0) / 0, and a2VF_(T,0) = 0. 

Hence, we must find r with Im r < such that 

r+oo 

/ X_ (r, s) ds / 0, and d^X^ (r, 0) = 0. 
Jo 

Moreover, we know from Step 3, section [3] that for Zmr < 0, the constraint ff^°°X 7^ is 
satisfied. Finally, the condition (SC) comes down to: 

dzX-. (r, 0) = 0, for some r with Im r < 0. 

To check this, and get a value for r, one can use a shooting method. For any r and any 
zq ^ 1, one can start from the approximation 

— IT 7 ]^ —IT 3 ]^ 

X_(r,zo) « 2g^p(__^^2)^ 9^X_(t,zo) ~ -AV^ ^exp(--Az^), 

and integrate backwards (jS.ip using a Runge-Kutta scheme. This gives access to the function 
dzX^{T, z), for any r and any z < zq. Then, a Newton- Raphson procedure allows to find a 
zero in {Imr < 0} for the function r dzX^^T, 0). Using such a procedure, we have found 

r « -0.706 - 0.706 i. 

Note that this value is proportional to 1 + i, as expected from the analysis. 

5.2 Simulation of the instability mechanism 

To observe the instability mechanism described in section [21 we have performed direct simu- 
lations of system ()2.ip . More precisely, we have considered the velocity 

Us{t,y) = Us(y) := 2y exp{-y^) 

(already studied in [B] in the inviscid case), and solutions of the type 

Ue{t,x,y) = ie''~'''dyVe{e-\y), v,{t,y) = e~^V,{e-H,y). 

The profiles = Vs{0,y) satisfy the singular perturbation problem 

idg + iUs)dyV, -iU% - ee''~'''dlV, = 

on Ve{9,y). One more differentiation gives the parabolic like equation 

{de + iUs) d^Ve - i Ui% - ed^yVe = 0, 

fulfilled with the boundary conditions 

Ve\y=0 = dyVe\y=Q = dyVe\y=0 = 0. 
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Figure 1: The correction to the inviscid eigenvalue. Plot of (^"""^(e) + Us{a)), seen as a 

function of the tangential frequency k = e~^. Top: real part. Bottom: imaginary part. As 
expected from the theory, both approximately converge to —0.9, as k goes to infinity. 




We have discretized this equation in space using finite differences on a stretched grid, and in 
time through a Crank-Nicholson scheme. Starting from initial random data {i.e. with random 
values at each gridpoint), we have computed its time evolution for values oi k = e"^ ranging 
from 1 to 3.10''. For sufficiently large times, one observes that the numerical solution y^^"^ 
behaves like 

in the sense that 

y ) ■ A6'y™"^(t,?/) 

gets independent of 6 and y. Computations show a relation of the type 

uj^^"^{e) Us{a) + Ve(-0.92-0.9H), 

see figure 1. Here a = is as usual the critical point of Us- This relation is in very good 
agreement with the theoretical prediction. 



UJ' 



■= -Us{a) + e^/^ ' ^ -Us{a) + ^/^ (-0.92 - 0.92 



It/"., fa) 1^/2 



r 



V2 

if we take for r the value —0.706 — 0.706i found in the previous subsection. Moreover, with this 
value of r, one can compute directly the solution V of the shear layer equation. After proper 
rescaling, this allows for comparison between the "numerical" and "theoretical" eigenmodes. 
More precisely, using the notations of (j2.5p . one can compare the functions 



..th ._ 1 ( v^'iy) My) \ _ 1 f vT'^iy) My) 



vZM ■■= -7= -^ear^ - and vZriy) = -7= 



^/E\vl^^{00) Va{oo)J V^r^loo) Wa(oo), 

which should both describe the correction to the inviscid eigenmode outside the shear layer. 
As regards the shear layer, one can compare 

vtiz) := ^"-^^^^^^^^^ and v^^iz) := ^ "^f + 
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Figure 2: Plots of Vg'l^ (dashed line) and v^^l^ (full line), seen as functions of y, at e = 10"''. 
Left and right figures correspond respectively to the real and imaginary parts. They match, 
as expected, outside the shear layer. 



Illustrations of these comparisons are given in figures 2 and 3. Again, we obtain an 
excellent agreement. This confirms that the instability mechanism we have described is indeed 
effective, and moreover dominates the linear dynamics ()1.6|) . 
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functions of z, at e = 10""^. Left and right figures correspond respectively to the real and 
imaginary parts. 
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Appendix : Well-posedness in the analytic setting 

We start from a simple estimate on the heat equation: 

For Uo e W"a'°°(M+), F G L^{0,T; Wa'°°{R+)), the solution U of 

dtU~dp = F on RxR+, U\t=o = Uq, U\y=o = 0, 
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satisfies 

^ ^ (ll^0||^y0,oo + 1 1 1 1^1(^/0,00-)) . (5 

This estimate follows directly from the representation formula 

U{t,y) = / S{t,y,z)Uo{z)dz + / / S{t - s,y, z)F{s, z) dz ds 
Jr+ Jo Jm.+ 

where the heat kernel in the half plane S{t, y, z) reads 

S{t, y, z) := G{t, y-z)- G{t, y + z), G{t, y) := exp(-2/V4t). 

Details are left to the reader. This estimate allows to prove proposition [TJ Indeed, 
decomposing 



the well-posedness is an easy consequence of the a priori estimate 

\\U''{t,-)\\w^,^ < Ce^'=*||C/'=(0,-)||vFO,c 
for some p. Now, the equation satisfied by is 

dtU^-dlU^ = ik (^U', j\''{t,z)dz - UsU^ 



Using (jS.ip . we get 



f/'(t,-)llwO- < C||C/'(0,-)llw0.oo + Csk / V'(^,-)llwO,oods 



where Cg depends on Ug- We conclude by the Gronwall lemma. 
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